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ABSTRACT 

We present VTL/ISAAC infrared spectroscopy of a sample of short orbital period 
cataclysmic variables which are candidates for harboring substellar companions. We 
have detected the K I and Na I absorption lines of the companion star in VY Aqr. The 
overall spectral distribution in this system is best fit with a M9.5 type dwarf spectra, 
implying a distance of 100 ± 10 pc. VY Aqr seems to fall far from the theoretical 
distribution of secondary star temperatures around the orbital period minimum. Fit- 
ting of the IR spectral energy distribution (SED) was performed by comparing the 
observed spectrum with late-type templates. The application of such a spectral fit- 
ting procedure suggests that the continuum shape in the 1.1-2.5 fim spectral region in 
short orbital period cataclysmic variables may be an useful indicator of the compan- 
ion spectral type. The SED fitting for RZ Leo and CU Vel suggests M5 type dwarf 
companions, and distances of 340 ±110 and 150 ± 50 pc, respectively. These systems 
may be placed in the upper evolution branch for short period cataclysmic variables. 

Key words: Stars: individual: VY Aqr, V485 Cen, RZ Leo, WZ Sge, CU Vel, 
HV Vir, 1RXSJ105010.3, Stars: binaries: close, Stars: Cataclysmic Variables, Dwarf 
novae, fundamental parameters, evolution, Stars: binaries: general, Stars: low-mass 



1 INTRODUCTION 

Cataclysmic variable stars (CVs) are semi detached bina- 
ries consisting of an accreting white dwarf and a red dwarf 
donor transferring matter to the compact object via the in- 
ner Lagrangian point. The orbiting gas interacts with itself 
dissipating energy by viscous forces, forming a luminous ac- 
cretion disk around the white dwarf. The spectral energy 
distribution of CVs in X-rays and ultraviolet is dominated 
by white dwarf and inner disk emission, whereas the accre- 
tion disk contribution is dominant in the optical and even- 
tually in the near infrared. In some cases the emission from 
the accretion disk in the IR is approximated by a power-law 
(Ciardi et al. 1998). However, in general we expect a much 
more complex emission from low-luminosity disks which may 
contain extended optically thin regions. On the other hand, 
the secondary star might contribute significantly in the in- 
frared. The determination of the secondary mass in CVs is 
key to understanding the secular evolution of these objects. 
Current theories state that the process of mass transfer be- 
comes linked with the loss of orbital momentum, so the bi- 

* Based on observations obtained at the European Southern 
Observatory, ESO proposal 65.H-0409(A). 
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nary period becomes shorter while the hydrogen secondary 
becomes less and less massive, eventually being eroded by 
the process, resulting in a kind of brown dwarf star when 
the orbital period approaches 80 minutes (e.g. Howell et 
al. 2001). An alternative scenario considers that most CVs 
may not yet have had time to evolve to their theoretical 
minimum orbital periods. In this case, for initial configura- 
tions with intermediate secondary masses, thermal-timescale 
mass transfer may occur. "If some nuclear evolution has al- 
ready occurred, this phase can shrink the binary drastically 
and strip the hydrogen-rich envelope from the donor, ulti- 
mately producing an ultrashort-period system with a low- 
mass, hydrogen poor and probably degenerate secondary" 
(from King & Schenker 2002). From the above it is evident 
that relevant observations for probing current theories of CV 
evolution should focus on the determination of the physical 
parameters of the secondary star for systems below the or- 
bital period gap. At present, four methods have been used to 
search for undermassive secondary stars in cataclysmic vari- 
able stars: (1) analysis of the spectral energy distribution us- 
ing multi- wave-band observations through the ultra-violet, 
optical, and infrared spectral regions (Ciardi et al. 1998; Ma- 
son 2001), (2) looking for signatures of the secondary star 
(Steeghs et al. 2001, Littlefair et al. 2000, Dhillon & Marsh 
1995), (3) "weighting" the secondary star in systems where 



© 1994 RAS 



2 R.E. Mennickent and M.P. Diaz 



the superhump and orbital periods are known (the primary 
mass is usually assumed, see for instance Patterson 2001), 
and (4) by spectroscopic diagnostic of the stellar masses in 
systems where the white dwarf is revealed by their optical 
absorption wings (Mennickent et al. 2001). Concerning the 
first two methods, one should mention that the determina- 
tion of basic properties of secondaries in CVs by comparison 
of their spectra with field stars calibrations is intrinsically 
uncertain. These results are prone to illumination and heat- 
ing of the companion photosphere. In addition, the absorp- 
tion spectrum may be affected by the filling of some lines 
with emission components. 

While the IR spectra of CVs has been measured and 
modeled in the past years there are only a few spectropho- 
tometric observations of low-luminosity systems in J,H and 
K band. In this paper we describe the spectra of 7 CVs with 
orbital period below the period gap. These candidates for 
systems at late evolutionary stages were selected for obser- 
vation using the following criteria: i. their short orbital pe- 
riod and/or ii. a low mass transfer rate, inferred from their 
high- amplitude dwarf-nova outbursts with long recurrence 
time. In the next section we describe the IR spectroscopic 
observations while in Section 3 a description of each spec- 
tra is given. A brief discussion of the observational results is 
made in Section 4. A few conclusions and perspectives are 
outlined in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

The infrared spectroscopic observations reported in this pa- 
per were obtained at ESO with VLT-Antu using ISAAC 
spectrograph in service mode. Data were taken under pho- 
tometric atmospheric conditions. An observing log is given 
in Table 1. Spectra in the J, H and K bands were obtained. 
Sufficient overlapping was assured for composing a single 
spectra ranging from 1.09 to 2.57 microns with a FWHM 
resolution of 12 (J) to 27 (K) angstroms. The data were 
reduced using IRAF0 by first applying the combined dark 
and flatfield images as supplied by the ESO service mode op- 
eration group. Median sky frames were then combined for 
each object and spectral window. This process made use of 
jittered images where the object is located at different po- 
sitions along the slit. Wavelength calibration was achieved 
by measuring the location of atmospheric OH emission lines 
(Rousselot et. al. 2000) in the sky background. The flux cali- 
bration was performed using observations of the AO standard 
HD216009, made with the same instrumental setup (but 
with a wider slit) by the ESO operation team as part of the 
service mode program. Finally, the telluric absorption fea- 
tures were corrected with the aid of the absorption template 
constructed by dividing the spectrum of HD216009 with a 
low order polynomial fit, excluding a few stellar absorption 
lines. The IRAF task "telluric" was used to find the best 
scale and shift factors which, when applied to the normal- 
ized telluric template, provided a reasonable correction of 

t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 



the telluric absorptions in HD216009 and science exposures. 
This procedure worked well, except for the regions between 
1.35-1.44 microns and 1.80-1.94 microns, characterized by 
heavy telluric absorption. These regions, corresponding to 
the ends of the J, H and K spectra, were excluded from the 
analysis and are not shown in this paper. Synthetic J, H, 
and K photometry of our calibrated standard star observa- 
tions were compared to broadband photometry by Carter 
& Meadows (1995) showing differences below 0.08 magni- 
tudes. Slit losses for our objects were estimated by assuming 
a Gaussian seeing profile centered in the slit aperture. For 
that we used the seeing value included in the report of the 
observing block and neglected any factor due to telescope 
guiding. These corrections were included in the J, H and K 
magnitudes and in the final spectrum as well. We are confi- 
dent that this method worked well, since our spectrophoto- 
metric magnitudes compare well with previously published 
photometry. In addition, we obtain a good match in the 
overlapping region between J,H and K spectra. 



2.1 Infrared spectral fitting 

An attempt to quantify the properties of the IR spectra 
was made by employing a numerical fitting procedure. The 
spectral energy distribution (SED) in the IR was tentatively 
parameterized by adding the contribution from a late-type 
template spectrum and power-law component. Although the 
emission of the disk should differ substantially from a power- 
law in the IR, we introduced this component as a first ap- 
proximation to the accretion disc continuum. 

While the power law is smooth in our wavelength do- 
main and basically affects the slope of the IR continuum, the 
detailed shape of the synthetic continuum in the J, H and K 
is strongly dependent on the stellar template contribution. 
A nonlinear least squares fitting procedure was calculated 
using the following equation: 

5(A) = a x T(A) + b x \ c (1) 

where S is the observed spectrum, T the red dwarf template 
spectrum, A the wavelength in microns and a,b,c parame- 
ters to be found. The parameters a, b and c were adjusted to 
minimize the reduced Chi-square between the observed spec- 
trum and a model fit. Of course, a and b are constrained to 
the positive domain. The data fit in this way were selected 
carefully avoiding emission lines and deep telluric bands. A 
sequence of template spectral types between Ml and L7 at 
2 subtypes step was taken from Leggett et al. (2001). The 
resolution of the template spectra (typically between 20—35 
A) was matched by convolving our data with a Gaussian 
with the appropiate FWHM. Due to the low velocity of the 
templates, it was not necessary to account for rotational 
broadening when fitting the spectrum. 



3 RESULTS 

The infrared spectra for all program stars are shown in 
Fig. 1. The synthetic magnitudes are given in Table 2 and 
the emission line parameters in Table 3. Table 4 summarizes 
the results of the spectral energy distribution fitting, giving 
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the parameters for the best fitting function as defined in 
Eq. f. 

3.1 VY Aqr 

VY Aqr is a cataclysmic variable showing one of the 
largest outburst amplitude among dwarf novae (Downes et 
al. 2001). Accordingly to current models for dwarf nova 
outburst (Osaki 1996), this is consistent with a very low 
mass transfer rate system. Spectroscopic studies in the op- 
tical region have revealed the orbital period (0'?06309(4) 
Thorstensen 1997), but not the stellar masses, likely due 
to the biasing nature of the emission line radial velocities 
(e.g. Augusteijn 1994). J-band spectroscopy by Littlefair et 
al. (2000) revealed spectral features of the secondary star 
but too weak to make an estimate of the spectral type. Ac- 
cording to the current CV evolution scenario (e.g. Howell 
et al 2001), the orbital period and the low mass transfer 
rate of VY Aqr should suggest a system beyond the orbital 
period minimum, probably containing an undermassive sec- 
ondary star. This is also supported by the difference between 
observed superhump period and the orbital period (Patter- 
son 2001). In this section we present direct evidence for this 
view, based on the unambiguous detection of the secondary 
star in the spectrum of VY Aqr. 

We measured synthetic magnitudes J = 15.17, H = 
14.66 and K = 14.14 with estimated errors of 0.10 mag, indi- 
cating that VY Aqr was observed in quiescence. Our J and K 
magnitudes compares well with the 15.24 ± 0.15 and 14.42 
± 0.13 values given by Sproats et al. (1996). The infrared 
spectrum of VY Aqr shows Bracket and Paschen emission 
lines. Fig. 2 clearly shows the K I doublets at 1.169-1.177 and 
1.244-1.253 and the Nal line at 1.141, which are signatures of 
a cool secondary star, confirming previous indications found 
by Littlefair et al. (2000). 

From the visual inspection of the spectral features de- 
tected in VY Aqr, we conclude that the secondary in VY 
Aqr is of spectral type later than M6 because of the follow- 
ing features: 1- the depth of K I doublets increases for later 
types. Even disregarding veiling, their depth suggests a type 
cooler than M5. 2- The depth of water band at 1.32 microns 
is too shallow in dwarf spectra earlier than M6. 3- Strong Al I 
lines at 1.313 and 1.32 fim, usually visible in dwarfs hotter 
than M6, are not seen in the object spectrum. On the other 
hand, a secondary spectrum earlier than L5 is suggested by: 
1- the strength of Nal 1.141 relative to the Kl doublet at 
1.169, 1.177 and 2- the depth of the water band. 

When fitting the spectral energy distribution (Fig 3), 
we observed that spectral types earlier than M7 fail to repro- 
duce the depth of K I lines in the J band and the continuum 
in the K band. On the other hand, spectral types later than 
L3 do not fit well the H and K band continuum shape. These 
cool types present a well defined CO bandhead at 2.29 )im 
which is not seen in our data. Our fit with a M9.5 type 
secondary is slightly better than that for M7 and L3 type 
templates, giving a \ 2 parameter about 15% lower. If such 
a spectral contribution is in fact due to the emission of the 
secondary star in the system one may estimate its temper- 
ature. Using the effective temperatures for L type dwarfs 
derived by Leggett et al. (2001) using structural models, we 
find T 2 = 2300 ± 100 K for the secondary star in VY Aqr. 
We note that our result is in apparent disagreement with 



the result of 1400 ± 50 K obtained by Mason (2001) by 
modeling the (non-simultaneous) UV-optical-infrared spec- 
tra. However, our result is still below the 2600 K upper limit 
she found by fitting her photometric dataset. Also our result 
is consistent with the upper limit suggested by Littlefair et 
al. (2000). The best fit with a M9.5 companion is shown for 
illustration in Fig. 3. Representative fits using types between 
M7 and L3 indicate that the secondary star may contribute 
with 55% to 45% of the flux at 2.17 /im, depending on the 
spectral type. Later types yield better fits for smaller flux 
fractions. Using the distance values for our templates from 
M7 to L2 (LHS3003 and LHS429 by van Altena et al. 1995 
and Kelu-1 from Dahn et al. 2000) and the flux fractions 
derived from the spectral fitting we were able to derive a 
distance estimate for the system between 80 and 120 pc, 
with a most likely value of 100 ± 10 pc. The spectroscopic 
parallax given above is in agreement with the distance of 
110 pc found by Augusteijn (1994) using the average abso- 
lute magnitude value for dwarf novae in outburst. 

3.2 V485 Cen 

This is a dwarf nova with an ultra-short orbital period of 59 
minutes (Augusteijn et al. 1996) . These authors found strong 
evidence for a significant contribution from the secondary to 
the spectrum at wavelengths longer than 5900 Angstroms, 
which was confirmed by Munari & Zwitter (1998), based 
on the observed color progression. According to Augusteijn 
et al. (1996), the most likely explanation for the extremely 
short orbital period, it is a non-fully degenerate secondary 
star with a low, but finite, hydrogen content. Considering 
our M-L spectral type templates, the best fit was reached 
with a simple power continuum with index -1.7. 

3.3 RZ Leo 

The long cycle length and large amplitude outburst relate 
this star to WZ Sge, however the orbital period (0.07651 
d, Mennickent & Tappert 2001) is longer and the spectral 
distribution redder (e.g. Mennickent et al. 1999). This fact, 
along with the rather normal mass ratio estimated from the 
superhump period (Ishioka et al. 2001) point to a normal 
secondary for this object. 

Our data seems to confirm this view. Fig. 4 shows that 
the best SED model for RZ Leo is reached with a M5 type 
spectrum and almost zero contribution from the power law. 
M3 and M7 subtype templates failed to match the shape of 
the spectrum, giving, in both cases, a \ 2 three times larger 
than the M5 template. We do not see the absorption lines of 
the secondary in the J-band spectrum probably due to our 
poor S/N in this range. Using Mj for M dwarfs from Leggett 
et al. (2000), we estimate a distance of 340 ±110 pc. Our 
J and K magnitudes, viz. 16.57 and 15.45, are similar to 
that given by Sproats et al. (1996), namely 16.56 ± 0.19 
and 15.55 ± 0.17. 

3.4 WZ Sge 

The spectral type of the secondary in this system has been 
constrained to be later than M7.5 (Littlefair et al. 2000) 
whereas an upper limit for the secondary star temperature 
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of 1700 K has been obtained by Ciardi et al. (1998). These 
spectral signatures along with the extremely small difference 
between superhump and orbital period (Patterson 2001), 
strongly suggests that the star is placed beyond the orbital 
period minimum. The discovery of emission lines of the sec- 
ondary during the 2001 superoutburst placed an upper limit 
on the secondary mass of 0.11 M©, but did not rule out 
the hypothesis of a non degenerate mass donor (Steeghs et 
al. 2001). 

Our infrared spectroscopy, obtained outside eclipse, 
shows a rich emission line spectrum, as previously found 
by Dhillon et al. (2000), Littlefair et al. (2000) and by Ma- 
son et al. (2000). We identify Paschen j3 double emission 
with a half-peak separation of 610 km/s and FWZI of 3900 
km/s. The Bracket series is also present as double emission, 
with the violet peak larger than the red one, except in the 
Br7 line where the asymmetry reverses. Our J magnitude, 
14.25, compares well with the J magnitude 14.2 ± 0.2 given 
by Ciardi et al. (1998). Our H magnitude, 14.07, is within 
the uncertainties of the 13.8 ± 0.2 value given by Ciardi 
et al. (1998). Our K magnitude, 13.78, is lower than the 
13.3 ± 0.2 value reported by Ciardi et al. (1998). This is 
the only case where a significant difference exists between 
one of our spectrophotometric magnitudes and previously 
published values. 

A simple power law provides the best fit for this object 
(Table 4). By looking at the scaled template line strength 
in the K band we estimate that the late type component 
contributes with less than 25 % to the flux in the K band 
for the range M9.5-L7. By assigning such upper limit to the 
companion star emission and using the astrometric parallax 
of 45 pc (J. Thorstensen 2001, private communication), we 
estimate Mk > 12 for the secondary. This figure corresponds 
to a secondary with a mass lower than or equal to the hy- 
drogen burning minimum mass (0.07 MQ), for ages between 
0.1 and 10 Gyr, accordingly to the models of very-low-mass 
stars and brown dwarfs with dusty atmospheres by Chabrier 
et al. (2000). 



3.5 CU Vel 

This object is a rather typical dwarf nova with orbital period 
of 0.0785d and moderate outburst credentials (e.g. Mennick- 
ent & Diaz 1996). The Lyman alpha absorption has been 
modeled as photospheric emission from a white dwarf with 
effective temperature of 18500 K (Gansicke & Koester 1999). 
Fig. 5 shows that the best SED model for CU Vel is reached 
with a M5 type spectrum and almost zero contribution from 
a power law continuum. As in the case of RZ Leo, M3 and 
M7 subtype templates failed to match the shape of the spec- 
trum, giving, in both cases, about three times larger than 
for the M5 template. We do not see the absorption lines of 
the secondary in the J-band spectrum probably due to our 
poor S/N in this range. The M5 spectral type is consis- 
tent with the 0.15 Mq secondary derived from the dynami- 
cal analysis of the Ha emission line by Mennickent & Diaz 
(1996). Following the same procedure as for RZ Leo, we find 
a distance of 150 ± 50 pc. 



3.6 HV Vir 

This is the faintest target in our sample, with magnitudes 
18.7 and 17.0 in J and K respectively. Our 30 min inte- 
gration time resulted in a very noisy spectrum. However, 
the spectrum seems to indicate line emission in Pa/3, and 
a rather flat continuum. The object has been linked to the 
sparsely populated group of WZ Sge stars (e.g. Leibowitz et 
al. 1994). Our SED algorithm reveals that a L5 a late type 
template contributing with about 39% to the total light at 
2.17/im gives a slightly better fit than a simple power law. 
With the addition of a L5 template we get a x 2 parame- 
ter slightly lower than for L3 and L7 templates. The \ 2 i s 
larger for power laws combined with other spectral types. 
The differences found in \2 are so small that we should 
quote the possible detection of a L5 type spectral contri- 
bution in HV Vir with extreme caution. Using the effective 
temperatures of L type dwarfs by Leggett et al. (2001), one 
find T eff = 1550 ± 200 K for the secondary of HV Vir, and 
a distance of 175 ± 15 pc. 

3.7 1RXSJ105010.3 

This ROSAT source was recently identified as a cataclysmic 
variable star, with an optical spectrum closely resembling 
WZ Sge, and likely containing a sub-stellar secondary (Men- 
nickent et al. 2001). The infrared spectrum is very similar to 
that of WZ Sge too, both in the shape of the continuum as 
well as in the rich emission line spectrum and line strengths. 
No traces of the secondary spectral features were found. 

As in the case of WZ Sge, our SED models reveal than a 
simple power law provides the best fit. This fact may suggest 
that, in spite of the low mass transfer rate, the accretion disk 
is more luminous than the secondary even in this spectral 
region. 



4 DISCUSSION 

In Fig. 6 we show the J-H:H-K color diagram for short- 
period dwarf novae, along with data for M-L type dwarfs. 
Colors for pure power law energy distribution are also in- 
cluded for comparison. We observe that the CV colors are 
distributed in a rather ordered way, following a well-defined 
track. Longer period systems (RZ Leo, CU Vel) show redder 
colors, above the region of the stellar photospheres, whereas 
systems around the orbital period minimum are bluer, and 
best represented by a simple power law continuum with in- 
dexes between -2.8 and -1.8. 

The most prominent emission lines are Paschen and 
Bracket lines. Helium lines were not observed in our 
(small) sample. Among the observed spectra, WZ Sge, 
1RXSJ105010.3, and VY Aqr show the strongest emission 
lines. These three systems are supposed to sustain very low 
mass transfer rates during quiescence, as may be inferred 
from their long recurrence times and large amplitude out- 
bursts. 

In this paragraph we compare our results on the sec- 
ondary star of short orbital period CVs with the evolution- 
ary predictions of the population synthesis code by Howell 
et al. (2001). At this point we emphasize the caveat already 
made concerning the reliability of secondary temperature 
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estimates from the observed emission of a particular spec- 
tral type. The Fig. 7 shows a comparison of the T e // — P or b 
CV evolution tracks near the orbital period minimum with 
our data and some additional data taken from the litera- 
ture. From this figure we conclude the following: 1- HV Vir, 
WZ Sge, EF Eri, WX Cet, LL And and SW UMa seem to 
be post-orbital period minimum systems. 2- It is difficult 
to reconciliate the positions VY Aqr in the diagram with 
the code's predictions. 3- In the same context, RZ Leo and 
CU Vel should be evolving toward the orbital period mini- 
mum. 

The results from the application of the spectral fitting 
procedure suggests that the infrared continuum shape in 
short-period cataclysmic variables may be an useful indica- 
tor of the companion spectral type. This point is specially 
important if we consider that, due to the limitation imposed 
by the spectrum S/N in such faint systems, it is not always 
possible to detect the secondary star individual lines, but 
do the continuum shape. Also, the method have the ad- 
vantage of avoiding the uncertainties associated with non- 
simultaneous multi-wavelength observations, although their 
predictive power might be not so good as the ideal case 
of modeling of simultaneous multi-wave-band observations 
(e.g. Mason 2001). In the future, we plan to observe a larger 
sample of CVs which are candidates for harboring brown 



dwarf 1: ke secondaries, in this regard, the application of the 
method to below-tne-gap polars at deep quiescence turns to 
be especially attractive (Howell & Ciardi, 2001). In these 
cases, the accretion of matter by the white dwarf virtually 
stops, resulting in the easier detection of the secondary star 
at infrared wavelengths. In addition, a more realistic descrip- 
tion of the disk emission in the IR will be incorporated to 
the analysis. 



5 CONCLUSIONS 

• We have found some evidence for a M5 type spectral 
distribution in RZ Leo and CU Vel. This may suggest that 
these systems are in the upper branch of the CV evolution 
track, i.e. prior to the orbital period minimum; 

• We find marginal evidence for a L5 type secondary star 
in HV Vir. A secondary star with spectral type M9.5 may be 
identified in VY Aqr. These systems have probably passed 
beyond the orbital period minimum. However, the position 
of VY Aqr in the T2-P rb diagram conflicts with the results 
of the population synthesis code by Howell et al. (2001). 

• For some objects, namely V485 Cen, WZ Sge and 
1RXSJ105010.3, we found no significant improvement of the 
spectral fitting by adding a stellar atmosphere template. 
This may indicate that even for such low luminosity sys- 
tems the accretion disk spectrum dominates the flux in the 
IR. 
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Table 1. Journal of Observations. The total integration time, in 
seconds, is given for every spectral band. 



Object 


UT-Date 


J 


H 


K 


VY Aqr 


23/05/00 


300 


600 


1200 


V485 Cen 


21/05/00 


600 




600 


RZ Leo 


21/05/00 


840 


840 


840 


WZ Sge 


20/05/00 


240 


480 


720 


CU Vel 


22/05/00 


120 


240 


480 


HV Vir 


19/05/00 


1800 




1800 


1RXSJ105010.3 


21/05/00 


240 


480 


720 



Table 3. Equivalent width (W(\) in A) and peak flux ( F(X) in 
units of 10 — 16 erg s — 1 cm -2 A -1 ) of the main emission lines. 



Object 


W(Pa^) 


F(P a/3 ) 


W(Br 7 ) 


F(Br 7 ) 


VY Aqr 


-52 


4.79 


-59 


1.45 


V485 Cen 


-18 


1.01 






RZ Leo 


-34 


1.01 


-37 


0.37 


WZ Sge 


-144 


17.3 


-142 


2.76 


CU Vel 


-41 


6.44 


-26 


2.30 


1RXSJ105010.3 


-107 


1.00 


-142 


0.32 



Table 2. Synthetic magnitudes and normalization factors used 
in Fig. 1, in units of 10 -17 erg cm -2 s _1 A . 



Object 


.1 


H 


K 


factor 


VY Aqr 


15.17 


14.66 


14.14 


29.7 


V485 Cen 


16.47 




15.30 


9.16 


RZ Leo 


16.57 


15.68 


15.45 


8.89 


WZ Sge 


14.25 


14.07 


13.78 


74.0 


CU Vel 


14.74 


13.94 


13.51 


46.6 


HV Vir 


18.68 




16.98 


1.06 


1RXSJ105010.3 


17.07 


16.77 


16.07 


5.20 



Table 4. Summary of the best SED models. We give the best 
template spectral type (Sp2), their relative contribution at A 2.17 
fim (in %) and the parameter c. N/A means not applicable. 



Object 


Sp2 


rel. cont. 


c 


VY Aqr 


M9.5 


44 


-2.2 


V485 Cen 


N/A 





-1.7 


RZ Leo 


M5 


100 


N/A 


WZ Sge 


N/A 





-2.7 


CU Vel 


M5 


100 


N/A 


HV Vir 


L5 


39 


-0.9 


1RXSJ105010.3 


N/A 





-2.1 



Figure Captions 

Fig.l Combined J,H,K spectrum for selected short or- 
bital period CVs. Fluxes are normalized by the factors 
shown in Table 2. 

Fig. 2 Comparison of the normalized spectrum of VY 
Aqr with the late-type templates. 

Fig. 3 The infrared spectrum of VY Aqr and the best 
composite SED fit (thick line). The individual template 
spectrum and power law continuum are shown. The spec- 
tra are scaled accordingly to the factor listed in Table 2. A 
zoom around the Kl doublets is also shown in the upper 
left graph. Main emission lines and absorption bands were 
excluded from the fit. 

Fig. 4 The infrared spectrum of RZ Leo and the best 
SED fit (thick line), corresponding to a M5 type template. 
The spectra are scaled to the flux at A = 1.7 /tm. 

Fig. 5 The infrared spectrum of CU Vel and the SED 
best fit (thick line), corresponding to a M5 type template. 
The spectra are scaled to the flux at A = 1.7 /an. 

Fig. 6 The J-H:H-K diagram for short-period CVs (solid 
triangles), M1-M6.5 dwarfs (open triangles, from Leggett 
et al. 2000) and late-M and L-dwarfs (open circles, from 
Leggett et al. 2001 and Reid et al. 2001). Data points corre- 
sponding to pure power law models are labeled accordingly 
to the power index and they are interpolated by a dotted 
straight line. CV data are from Table 2. 

Fig.7 Cataclysmic variables close to the orbital period 
minimum. Observations are compared with results of the CV 
population synthesis code by Howell et al. (2001). We have 
used spectral type - temperature calibrations based on data 
of M-L dwarfs by Leggett et. al. (2000, 2001). Data for LL 
And and EF Eri are from Howell & Ciardi (2001), for WX 
Cet, EF Peg and SW UMa from Mason (2001) and for WZ 
Sge (a temperature upper limit) from Ciardi et al. (1998). 
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